Photoproduction of Quarkonium in Proton-Proton and Nucleus-Nucleus Collisions 
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We discuss the photoproduction of T and J/ip at high energy pp, pp and heavy ion coUiders. We 
predict large rates in pp interactions at the Fermilab Tevatron and in pp and heavy-ion interactions 
at the CERN LHC The J/tp is also produced copiously at RHIC. These reactions can be used 
to study the gluon distribution in protons and heavy nuclei. We also show that the different CP 
symmetries of the initial states lead to large differences in the transverse momentum spectra of 
mesons produced in pp vs. pp collisions. 
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Photoproduction has traditionaUy been studied with 
fixed target photon beams, at electron-proton colliders, 
and, to a limited extent, at relativistic heavy ion collid- 
ers. However, energetic protons also have large electro- 
magnetic fields, and high energy pp and pp colliders can 
be used to study photoproduction, at photon energies 
higher than are currently accessible. These photopro- 
duction reactions are of interest as a way to measure the 
gluon distribution in protons at low Feynman— a:. 

In this letter, we study photoproduction of heavy quark 
vector mesons in pp collisions at the Fermilab Tevatron 
and in pp collisions at the Relativistic Heavy Ion Col- 
lider at Brookhaven and at the Large Hadron Collider 
(LHC) at CERN. At the Tevatron, the high rates allow 
for detailed measurements of gluon distributions around 
a; « 1.5 - 5 X IQ-^, and at the LHC x^2-7x lO"" can 
be reached. We also discuss the coherent photonuclear 
production of T in nucleus-nucleus collisions at RHIC 
and the LHC. 

We show that the different CP symmetry between 
the pp and pp initial states leads to large differences in 
the transverse momentum, pT spectra of the produced 
mesons. Finally, we discuss how these events can be sep- 
arated from purely hadronic interactions. 

We use data from HERA and fixed target experiments 
on exclusive photoproduction of heavy vector mesons in 
photon-proton interactions [H as input to our calcula- 
tions. Because data on the T is limited, we use QCD 
based models as a basis for parameterizations of the cross 
sections. The paucity of experimental data on photopro- 
duction of the T leads to a relatively large uncertainty 
in the parameterized cross sections, but is also a strong 
motivation for investigating new production channels. 

The total J/^* photon-proton cross section for quasi- 
real {Q^ ~ 0) photons has been measured from near 
threshold up to photon-proton center-of-mass energies, 
Wjp, up to 200 GeV. The cross section increases with 
Wjp roughly as Wi^^. We parameterize a^p{W^p) = 
1.5 • W^-^ [nb] (W in GeV). A drawback of this parame- 



terization is that there is a discontinuity at the threshold 
energy, Wjp = rrip + nij/^. 

The two measurements of the T at HERA both have 
significant uncertainties. The Zeus collaboration mea- 
sured a ■ Br(T fifi) — 13.3 ± 6.0(stat.)l2'3(syst.) pb 
at a mean center-of-mass energy of {W^p) —120 GeV 
0. The HI collaboration found a ■ Br{T iifj,) = 
19.2 ± 9.9(stat.)±4.8(syst.) pb at {W^p) =143 GeV 0. 
Both experiments estimate that roughly 70% of the sig- 
nal comes from the T(IS') state. 

The leading-order expression for the photoproduction 
of a vector meson of mass My is (a 



da{'^p Vp) 



dt 



3aM, 



V 



-167T^xg{x,M^/4f . (1) 



Two more sophisticated calculations have considered the 
use of relativistic wave functions, off-diagonal parton dis- 
tributions, and NLO contributions Although the 
approaches differ, the final results are in good agreement. 
The cross section for T(1S) production scales roughly 
as W^J. We use a parameterization which is consistent 
with both HERA results: a^piW^p) = OM-W^J [pb] (W 
in GeV) . An alternative method based on parton- hadron 
duality gives cross sections ^30-50% larger depending on 
Wjp Q. Our calculations are for the T(1S). 

We estimate the uncertainties in the T cross section by 
fitting the HI and Zeus data to the function cF-yp(W-yp) — 
C-W^p. The constant C is determined from two fits, one 
with the experimental errors (quadratic sum of statistical 
and systematical) added to the measured value and the 
other with the experimental errors subtracted. These fits 
give C = 0.175 [pb] and C = 0.054 [pb], respectively. 

The cross section to produce a vector meson in a 
proton-proton collision is 



(t{p + p p + p + V) 



dn 
dk 



,{k)dk. (2) 



As will be discussed below, quantum mechanical inter- 
ference alters the transverse momentum (j)t) spectrum, 
but does not affect the total cross section significantly. 
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We use the photon spectrum from Ref . f?^ : 



RHIC p+p ^ p+p+J/f 



Tevatron p+p -» p+p+JW 



dn 
dk 



a 
27rfc 



1 . 11 3 



(3) 



where A = 1 + (0.71Gel/2)/g2^^^ and g2^,„ « (fc/7)^ It 
is derived using a proton form factor, F{Q^) — 1/(1 + 
(5V(0-71Gey2))2_ This spectrum is close to that of 
a point charge with a minimum impact parameter of 
bmin = 0.7 fm. It corresponds to emission of a photon 
with the proton remaining intact. For exclusive vector 
meson production, the target proton must also remain 
intact. This stricter requirement slightly decreases the 
effective photon flux. To estimate this uncertainty, we 
have also performed calculations with a photon spectrum 
corresponding to a point charge with a minimum impact 
parameter of &mm = 1-0 fm. 

The rapidity, y, of a produced state with mass My is 
related to the photon energy through y — \n{2k/Mv)- 
Using this relation in Eq. |21 and differentiating gives 



da dn 



(4) 



Interchanging the photon emitter and target corresponds 
to a reflection around y = 0; the total cross section is the 
sum of the two possibilities. The rapidity distributions 
are shown in Fig. The calculations are for collision 
energies of Vs = 500 GeV at RHIC, = 1-96 TeV at 
the Tevatron, and ^ U TeV at the LHC. The sohd 
and dashed histograms are for the parameterizations 
a^piW^p) = 1.5 • W°p^ [nb] and a^p{W^p) = 0.06 • W^J 
[pb] for the J/'^ and T, respectively. The solid his- 
togram is for the photon spectrum in Eq. 13 while the 
dashed histogram is for a point charge with a cut-off at 
b,nin = 1.0 fm. The grey band shown for the T corre- 
sponds to the two fits to the data described above (with 
the photon spectrum of Eq. The sharp cut-off at 

large rapidities for the J/^P is due to the discontinuity at 
threshold in the parameterization of a-yp. Using the pho- 
ton spectrum in Eq. |31 the total cross sections for the T 
are 12 pb, 120 pb, and 3.5 nb at RHIC, the Tevatron, and 
the LHC, respectively. For the J/tp, the corresponding 
cross sections are 7.0 nb, 23 nb, and 120 nb. 

Exclusive J/^ production in pp interactions was pre- 
viously considered by Khoze et al.^. They use a very 
different approach, based on the proton energy lost. An- 
alytical or numerical comparisons between their results 
and ours are therefore difficult. 

At the design luminosity for pp at the Tevatron (2 x 
lO^^cm^^g-i) the T production rate is 76/hr. Even 
with the 2.5% per flavor branching ratios to good 
signals should be observable at the Tevatron. For exam- 
ple, in 1 year (8,000 hr) at 33% overall efficiency, 5,000 de- 
cays should be observable in each of the e+e" and 
channels. At RHIC, the lower energy, luminosity and 
shorter pp running time decrease the signal. However, 




FIG. 1: Rapidity distributions for photoproduction of J/4' 
and T(1S) mesons in pp and pp interactions at RHIC, the 
Tevatron, and the LHC. The curves are explained in the text. 



with the planned RHIC II luminosity upgrade, T photo- 
production could be studied. The situation is better for 
the J/^*, where the production rate at RHIC is 250/hr. 
At the LHC, the J/^ and T rates are both high. 

The mid-rapidity photoproduction cross section for T 
at the Tevatron is about 0.1% of the hadronic inclusive T 
cross section |lOj . Similarly, the photoproduction cross 
section at mid-rapidity for the J/^ at RHIC is about 
0.1% of the hadronic inclusive cross section measured at 
^/s — 200 GeV (the total cross section is expected to be 
about twice as large at -^s = 500 GeV) [Tll| . 

Although the photoproduction cross section is a small 
fraction of the hadronic cross section, separation of this 
reaction channel seems possible given the very different 
character of the photon induced events. We will discuss 
some of the selection criteria, and estimate their effec- 
tiveness at rejecting hadronic events. 

Hadronically produced vector mesons have px My- 
In contrast, almost all of the photoproduced mesons have 
Pt < i GeV/c (cf. Fig. 2). A px < I GeV/c cut elimi- 
nates about 94% of the hadroproduced T at the Tevatron 
jlOj l , while retaining almost all of the photoproduction. 

As long as both protons remain intact, the vector me- 
son will not be accompanied by any other particles in the 
same event. In contrast, in hadronic events, the produced 
particles are distributed over the available phase space. 
If the average charged particle multiplicity is {duch/dy), 
then the probability of having a charged particle free ra- 
pidity gap with width Ay is exp(— Ay • {drich / dy)) . The 
mean particle densities at mid-rapidity are (drich/dy) = 
3.0 at = 500 GeV, (driah/dy) = 4.0 at = 1.96 TeV, 
and {drich /dy) ~ 5.5 at -^s = 14 TeV (extrapolated), ne- 
glecting any possible difference between p and p (12] . Re- 
quiring that the vector meson be surrounded by particle 
free regions (rapidity gaps) with a total width Ay = 3.0 
will reduce the background by a factor of « lO"** at 
RHIC, « 10-5 at the Tevatron, and w 10"^ at the LHC. 
The total width can be split between two or more gaps. 
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RHIC p+p ^ p+p+T Tevatron p+p ^ p+p+T 
60 



it corresponds to a charge-parity (CP) operation. Since 
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FIG. 2: da/dydt for photoproduced T(1S) mesons at mid- 
rapidity in pp and pp collisions at RHIC and the Tevatron. 
The inset has an expanded t scale. 



e.g. two gaps of width 1.5, or one of width 3.0. These 
gaps fit within the acceptance of existing and planned de- 
tectors, and provide more rejection power than is needed. 
The selection could be improved by using calorimetry to 
detect neutral particles in the gaps. 

The CDF collaboration has identified a sample of ex- 
clusive J ftps they do not give a cross section, but 
most of the yield is in the region p^j- < 1 GeV/c, as ex- 
pected for photoproduction. 

Exclusive vector meson production in pp and pp colli- 
sions differs from production in ep or eA collisions in that 
both projectiles can act either as target or photon emit- 
ter. For very small momenta of the produced state, one 
cannot distinguish which proton (or anti-proton) emitted 
the photon and which acted as target, so adding the cross 
sections is not justified. This interference was studied for 
nucleus- nucleus collisions 

The interference is best understood in a plane perpen- 
dicular to the direction of the beams. The impact param- 
eters for photoproduction are typically a few fm because 
of the long range of the electromagnetic force, but the 
production is always localized to one of the two projec- 
tiles. Because of symmetry, the effect is largest for y = 0. 

The differential cross section, da/dydt, where t is the 
momentum transfer from the target (t « pf^), may be 
written as an integral over the impact parameter, b: 



da 
dydt 



b>b„ 



\Ai + A2\U%. 



(5) 



Ai and A2 are the amplitudes for production off each of 
the two targets. At mid-rapidity, considered here, \Ai \ — 
IA2I. The amplitude is normalized to the cross section 
for a single source (Eq. 0}. 

If the produced vector mesons are treated as plane 
waves, Ai = Aq expiip ■ x), the total amplitude is 



1^1 



A.\'^2 



1^1 P (l±cos(pT-6) 



(6) 



The sign of the cos-term depends on the symmetry of 
the system. In a pp collision, moving the vector meson 
emission (scattering) from one proton to the other corre- 
sponds to a parity transform. In a pp collision, however. 



the vector meson has quantum numbers J 



PC _ 1 — 



1 — , the 



interference is destructive for pp and constructive for pp. 
The sign in Eq. El is "— " in pp collisions (as in nucleus- 
nucleus collisions [lj|) and in collisions. With ad- 
equate statistics, the interference might be used to search 
for CP violation. 

This interference alters the vector meson px spectrum 
near mid-rapidity. Without interference, the pr spec- 
trum is that for production off a single (anti-)proton. 
This spectrum is the convolution of the photon trans- 
verse momentum spectrum with the spectrum of trans- 
verse momentum transfers from the target For 
Pt > h/ {b), the cos{pt -byteim in Eq.^loscillates rapidly 
as b varies, and the net contribution to the integral 
will be zero. For small transverse momenta, however, 
Pt <C 1/(6), p ■ ^ ~ for all relevant impact parameters, 
and interference alters the spectrum. This is illustrated 
in Fig. 121 which compares da/dydt with and without in- 
terference at RHIC and the Tevatron. For Fig.|21 we im- 
posed a cut bmin — 1-0 fm; this has a small effect on the 
spectrum. The interference is large for t < 0.05 GeV^/c^. 
The different sign of the interference in pp and pp is 
clearly visible. 

In addition to pp and pp interactions, vector mesons are 
produced in coherent ultra-peripheral nucleus- nucleus 
collisions The STAR collaboration has observed 

Au + Au ^ Au + Au + p° at RHIC 16\. With a cut 
on Pt < lOOMeV/c, the signals were quite clean. 

For coherent production, the momentum transfer from 
the nucleus to the vector meson is determined by the 
nuclear form-factor. Since the T(IS') has a small cross 
section to interact with a nucleon, hadronic shadowing 
should be negligible for it. The forward scattering am- 
plitude scales with the number of nucleons. A, squared: 



da{-fA -> T^) 



dt 



A' 



da{'-fp — > Tp) 



dt 



t=0 



(7) 

A Woods-Saxon distribution is used for the nuclear form 
factor F{t). The total photonuclear cross section is 
the integral of Eq. [7| over all momentum transfers, 
t > tmin = [M|-/4fc7]^. The input to the calculation is 
again the parameterizations of the total photon-proton T 
cross section discussed above. When determining the for- 
ward scattering amplitude from the total photon-proton 
cross section, an exponential t-dependence is assumed 
with the same slope, 4 GeV~^, as for J/^* production. 
This leads to a forward scattering amplitude about 5% 
lower than if the proton form factor above had been used. 

With this, and the photon spectrum in TB'l, the to- 
tal cross section and the rapidity distributions can be 
calculated. Figure |31 shows da/dy for T production in 
Si-I-Si interactions at RHIC (7 = 135) and Pb-f Pb inter- 
actions at the LHC (7 = 2940). The total cross section 
(sohd curve) is 0.72 nb for Si+Si at RHIC and 170 ^ih for 
Pb+Pb at the LHC. 
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RHIC Si+Si ^ Si+Si+T LHC Pb+Pb ^ Pb+Pb+T 




FIG. 3: Rapidity distributions of T mesons produced in co- 
herent photonuclear interactions at RHIC and the LHC. The 
solid curves correspond to the parameterization a^p{W-,p) — 
0.06 ■ W^''' [pb], and the grey bands show the uncertainty in 

Because of the low photon flux, the cross section at 
RHIC is rather low. At design luminosity for Si+Si 
(4.4 • 10^^cm~^s~^), about 300 Ts are produced in a 
RHIC year (10^ s). The situation is better with Pb-ions 
at the LHC. The design luminosity (1 • lO^^cm^^s^^) 
corresponds to a production rate of about 0.02 Hz or 
roughly 60 Ts per hour. The experimental identifica- 
tion in heavy-ion interactions is relatively easy because 
of the coherence requirement, which limits production to 

PT < V^h/R [3. 

Eq. ^ shows that the forward scattering amplitude 
for T production is proportional to the gluon density 



squared. A 30% reduction in the nuclear gluon density 
would roughly halve the cross section. Photonuclear T 
production at the LHC should be a sensitive probe of 
nuclear gluon shadowing in the range x « 2 x 10"'^. 

The coherent production of T at the LHC was studied 
recently by Frankfurt et al. [iJl- Our result (solid curve 
in Fig. O is about 10% higher than their result for the 
impulse approximation (no shadowing). The difference 
may be due to the slightly different photon spectrum and 
slope of da/dt in photon-proton interactions. With nu- 
clear gluon shadowing, the cross section at mid-rapidity 
may be reduced by as much as 50% for Pb-|-Pb \V^ . 

Photoproduction of other final states should also be ac- 
cessible at existing and future pp and pp colliders. For ex- 
ample, photoproduction of open charm and bottom could 
be used to measure gluon distributions. These events 
would have only a single rapidity gap, but the experi- 
mental techniques should be similar. 

To summarize, we have calculated the cross sections 
for photoproduction of heavy vector mesons in pp and pp 
collisions. The cross sections are large enough for this re- 
action channel to be observed experimentally. The da/dt 
is distinctly different in pp and pp collisions because of the 
interference between the production sources. The cross 
section for producing T mesons in coherent photonuclear 
Pb-|-Pb interactions at the LHC is large. Because of the 
distinctive experimental signature, these reactions should 
be easy to detect. 

This work was supported by the U.S. Department of 
Energy under Contract No. DE-AC-03076SF00098. 
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